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Abstract. Recently, high resolution observations with the help of the near- infrared adaptive 
optics integral field spectrograph SINFONI at the VLT proved the existence of massive and 
young nuclear star clusters in the centres of a sample of Seyfert galaxies. With the help of 
high resolution hydrodynamical simulations with the PLUTO-code, we follow the evolution of 
such clusters, especially focusing on mass and energy feedback from young stars. This leads 
to a filamentary inflow of gas on large scales (tens of parsec), whereas a turbulent and very 
dense disc builds up on the parsec scale. Here, we concentrate on the long-term evolution of the 
nuclear disc in NGC 1068 with the help of an effective viscous disc model, using the mass input 
from the large scale simulations and accounting for star formation in the disc. This two-stage 
modelling enables us to connect the tens of parsec scale region (observable with SINFONI) with 
the parsec scale environment (MIDI observations). At the current age of the nuclear star cluster, 
our simulations predict disc sizes of the order of 0.8 to 0.9 pc, gas masses of 10 6 Mq and mass 
transfer rates through the inner boundary of 0.025 Mq /yr in good agreement with values derived 
from observations. 
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1. Introduction 

Nuclear activity is an important phase in the evolution of galaxies. Whenever enough 
gas is fed onto the central supermassive black holes, galaxies become active and their 
nuclear region lights up significantly, compared to non-active phases of their evolution. 
The infalling material builds up a fast rotating accretion disc, which heats up viscously. 
The emerging UV/optical radiation illuminates a gas and dust reservoir - the so-called 
molecular torus. Given the high opacity of the dust in this wavelength range, the torus 
morphology is able to geometrically unify two observed classes of active galaxies: type 1 
objects, where the torus is viewed face-on (and all of the characteristics of the central 
region are visible) and type 2 objects (torus viewed edge-on), where most of the light from 
the accretion disc is absorbed and reemitted in the infrared regime. This is the essence 
of the Unified Scheme of Active Galactic Nuclei. Despite its general success when probed 
with observational data, the nature of the obscuring gas and dust reservoir - the torus, as 
well as the triggering mechanism of such phases are still a matter of active research. To 
assess the final stages of mass transport through the central tens of parsec region, nearby 
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objects are needed, due to the limitations on resolution posed by currently available 
instruments. Therefore, Seyfert galaxies are an ideal testbed to study these processes. 
MIDI (MID-infrarcd Interferometer) observations recently resolved the innermost parts of 
the torus directly for the first time (e. g. Jaffe et aL(2004)||Tristram et aL(2007)| Tristram 



et aZ.(2009) Burtscher et oZ.(2009) I and found evidence for a clumpy dust distribution. 



In these proceedings, we numerically model a scenario, which is able to describe the build- 
up and evolution of a nuclear gas disc or torus. Our simulations are put into context in 
Davies et al. (these proceedings). 



2. Three-dimensional hydrodynamical simulations 

Recently, high resolution observations with the help of the near-infrared adaptive op- 
tics integral field spectrograph SINFONI at the VLT proved the existence of massive 



and young nuclear star clusters in the centres of a sample of Seyfert galaxies (Davies 



et al. (2007) and Davies et al., these proceedings). With the help of three-dimensional 



high resolution hydrodynamical simulations with the Pluto code, we follow the evo- 
lution of such clusters. The gas ejection of their stars provide both, material for the 
obscuration within the mentioned Unified Scheme of AGN and a reservoir to fuel the 
central, active region and it additionally drives turbulence in the interstellar medium 
on tens of parsec scales. We start our simulations after the very violent phase of super- 
nova type II explosions and fast winds. In order to enable a direct data comparison, we 
constrain our input parameters with observed values of the well-studied nearby Seyfert 
galaxy NGC 1068. A summary of the parameters used for the simulations shown in these 
proceedings is given in Table[l] After approximately 50Myrs following the starburst, 



Table 1. Parameters of our standard models. 



Hydrodynamical model Effective disc model 

Parameter Value Reference Parameter Value 
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Mass of the black hole (Mbh), normalisation constant of the stellar potential (M*), initial gas 
mass (Mg°a), cluster core radius (R c ), torus radius (Rt), inner radius (-Rin), outer radius (R ou t), 
stellar velocity dispersion (<r*), exponent of the angular momentum distribution of the stars (/?), 
initial gas temperature (Tini), normalised mass injection rate (M n ), mass of a single injection 
(Mpn) and adiabatic exponent (F), thickness of the disc (<5 = disc scale height / radius of the 
disc), gas temperature (T), alpha viscosity parameter (a v ), age of the nuclear star cluster at the 
beginning of the simulations (t^^ter) an d at the end (t^w. Br) an< ^ resolution of the simulation 



(fir 



The references are: L0 3|Lodato fc Bertin (2003) 
(Jungwiert et al. (200V 



D07 



fc Matthews(2003) and J01 



Davies et a 7 . (2007) G03 Gallimore 
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Figure 1. Snapshots of the density distribution in a meridional plane of our 3D hydrodynamical 
standard model after (a) 0.2 orbits (corresponding to 7- 10 4 yr) and (b) 1.7 orbits (approximately 
6 ■ 10 yr). Panel (c) shows the temperature distribution in the same meridional plane after 1.7 
orbits. Panel (d) is a zoom into the innermost part, showing the nuclear disc component. 



the mass loss of the newly born stellar population is dominated by slow winds and the 
ejection of planetary nebulae respectively. Then, the input of energy into the ambient 
medium is low enough to enable transport of gas towards the centre allowing for the 
triggering of activity there (see also Davies et al., these proceedings). A similar time 
delay between the nuclear starburst and the onset of nuclear activity has also been found 



observationally, leading to the same conclusions as our work (Davies et aZ.(2007)l. We 
model the mass-loss of the stellar population as the ejection of single expanding clumps 
of gas, which we give the velocity of the underlying stellar distribution at the point of 
emission, made up of a rotation component and a random component]!] These blobs of 
gas - confined by cooling instability - merge to form larger entities, thereby dissipating a 
fraction of their initial turbulent motions and get transported towards the centre, form- 
ing slightly elongated filaments. These processes lead to a vertically wide distributed 
clumpy or filamentary inflow of gas on tens of parsec scale (see Fig. ]l^i,b,c), whereas 
a turbulent and very dense disc builds up on the parsec scale (sec Fig. ]T]i). A similar 
two-component structure has been found in intcrfcromctric observations with the help 



of the MIDI instrument as well in NGC 1068 Raban et a?. (2009) and the Circinus galaxy 
Tristram et aZ.(2007)| 

Being computationally very expensive simulations and as we currently are unable to 
take all physical processes happening in the innermost few parsec into account in the 
three-dimensional hydrodynamical calculations, we treat the inner disc separately with 
the help of an effective disc model, as described in Sect. [3] 



3. Effective disc modelling 

In order to capture the relevant physics in the inner region, we treat this disc separately 
in a one-dimensional, axisymmetric simulation, which also enables a direct comparison 
with observations. Due to the generally unknown physical mechanism to generate an- 
gularm momentum and mass transfer in the disc, we concentrate on the effects of a 
parametrised turbulent viscosity and additionally take star formation into account. Fol- 

f These nuclear star clusters possess a significant velocity dispersion and rotate with sub-Ke- 
plerian velocity. 



122 



M. Schartmann et al. 



lowing |Pringle (1981)| and |Lin fc Pringle (1987)[ the viscous evolution of such a disc can 
be described by the following differential equation: 



dt 



T,(t,R) 



I A 
R dR 



dR 



(y a Yi(t,R) R 3d ^B\ 



dR J 



dR 



(R 2 n) 



= S, 



(3.1) 



where £ is the gas surface density within the disc, R is the radial distance from the 
centre, t the time, v a the assumed alpha viscosity, f2 the Keplerian rotation frequency and 
S includes a source term for the mass input from our three-dimensional hydrodynamical 
simulations (Sect. [2]) and a sink term, taking star formation into account with the help 
of the Kennicutt-Schmitt law. The infalling material is placed according to its angular 
momentum in a Keplerian rotating disc structure (Fig. |2^). Gas consumed into stars is 
removed from the simulations and will no longer participate in the dynamical evolution. 
We solve Equ. |3.l| numerically with the help of MATLAB'sjfJpdepe solver. This two-stage 
modelling process enables us to directly compare our results to observed properties. 



4. Results and data comparison 

Fig. [2j) shows the comparison of the mass contributions of the various components 
from the ID modelling for an assumed alpha viscosity parameter of a — 0.05 and a gas 
temperature of 400 K. The input parameters of the simulations have been constrained 
by observations of the nearby Seyfert 2 galaxy NGC 1068 (see Table[lJ. At the current 
age of the starburst in NGC 1068 of approximately 250 Myr, the largest fraction of 
the integrated gas mass introduced into the disc simulations (given by the dashed line in 
Fig. [2]d) has been accreted with the help of viscous processes (dash-dotted line) . A slightly 
smaller amount has been transformed into stars and a few times 10 6 M Q of gas remains 
in the disc (solid line), which has been built up within a few million years. The latter 

f http:/ /www. mathworks.com/products/matlab/ 




R [pc] cluster age [Myr] 

Figure 2. (a) Mass input into our effective disc model, shown as the growth of the gas sur- 
face density at different radii, flowing in through a sphere with radius of 2.5 pc from our 3D 
Pluto turbulent torus simulation. The input radii are set according to an angular momentum 
distribution of a Keplerian disc. Shown is the logarithm of a histogram for 100 timesteps with 
a resolution of approximately 3400 yr between an evolution time of 0.7 Myr and 1.0 Myr, where 
the 3D hydro simulation is already in an equilibrium state. The triangles denote the maxima 
of each radial bin, the parabola is the parametrisation used as source term in the effective disc 
simulations, (b) Comparison of mass contributions of the various components of our ID effective 
disc model for an alpha viscosity parameter of 0.05 and a gas temperature of 400 K. The vertical 
dashed line denotes the estimated current age of the nuclear star cluster in NGC 1068. 
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is in good comparison with the disc mass Kumar (1999)| finds by fitting a clumpy disc 
model to the maser detections. Fig. [3^ shows the surface density distribution for three 
simulations with varying strength of the viscosity: a — 0.05, 0.1 and 0.2. For comparison, 
the range of slopes of the density distribution as inferred from MIDI observations of a 
sample of nearby Seyfert galaxies is overplotted as dashed line segments. The sizes of the 
corresponding dust distributions we derive from these curves amount to 0.8 to 0.9 pc, 
which is remarkably similar to values obtained with the help of the MIDI instrument 
of approximately 0.7pc and to the estimated extent of the maser disc (0.65 to 1.1 pc). 
Finally, the total accretion rate through the inner boundary of our disc simulations is 
plotted in Fig. [3]d for the a parameter study discussed above. Larger values of a clearly 
lead to a larger mass accretion rate and at the current age of the nuclear stellar cluster, 
the derived values agree w ell with typic al accretion rates of Seyfert galaxies of a few 
times 10~ 3 to 10~ 2 M /yr ( |Jogee (2006)) . However, NGC 1068 seems to be in a heavily 
accreting state, which can be accomodated in our model only, when assuming e. g. clumpy 
accretion. 



5. Concluding remarks 

In these proceedings, we show that evolving stars from a massive and young nuclear 
star cluster, as found in nearby Seyfert galaxies provide enough gas to assemble a parsec- 
sized nuclear gas disc. To this end, we combine a three-dimensional treatment of the 
mass-loss from an evolving nuclear star cluster with a simplified model for the innermost 
parsec scale region, where a nuclear disc builds up. As far as possible, we derive input 
parameters of our simulations from observations of the nearby and well-studied Seyfert 2 
galaxy NGC 1068. This two-stage analysis enables us to (i) do a long term evolution study, 

(ii) link the tens of parsec scale region of galactic nuclei (observed with the SINFONI 
instrument) to the sub-parsec scales (probed by MIDI and in water maser emission) and 

(iii) test our model directly with a large number of observational results. At the current 
age of its nuclear starburst of 250 Myr, our simulations yield disc sizes of the order of 0.8 
to 0.9 pc, gas masses of 10 6 M Q and mass transfer rates of 0.025 M@/yr through the inner 
rim of the disc in good comparison with observed disc and torus properties. On basis 
of these comparisons, we conclude that the proposed scenario seems to be a reasonable 
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Figure 3. (a) Final gas density distribution at a cluster age of 250 Myr. Superposed as dashed 
lines are the two limiting curve s for the range of possib le surface density distributions as deter- 
mined from MIDI observations Kishimoto et ai!. (2009) (b) Total accretion rate of gas through 
the inner boundary of the domain. 
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model and shows that nuclear star formation activity and subsequent AGN activity are 
intimately related. 
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